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We present an improved many-body T-matrix theory for partially Bose-Einstein condensed 
atomic gases by treating the phase fluctuations exactly. The resulting mean-fleld theory is valid 
in arbitrary dimensions and able to describe the low-temperature crossover between three, two and 
one-dimensional Bose gases. When applied to a degenerate two-dimensional atomic hydrogen gas, 
we obtain a reduction of the three-body recombination rate which compares favorably with ex- 
periment. Supplementing the mean-field theory with a renormalization-group approach to treat 
the critical fluctuations, we also incorporate into the theory the Kosterlitz-Thouless transition that 
occurs in a homogeneous Bose gas in two dimensions. In particular, we calculate the critical con- 
ditions for the Kosterlitz-Thouless phase transition as a function of the microscopic parameters of 
the theory. The proposed theory is further applied to a trapped one-dimensional Bose gas, where 
we flnd good agreement with exact numerical results obtained by solving a nonlinear Langevin fleld 
equation. 



PACS numbers: 03.75.Fi, 67.40.-w, 32.80.Pj 

I. INTRODUCTION 

Low-dimensional Bose gases have recently attracted at- 
tention both experimentally and theoretically. The inter- 
est in these systems stems from the fact that the physics 
of low-dimensional systems is fundamentally different 
from the physics of systems in three dimensions. One 
and two-dimensional Bose-Einstein condensates have re- 
cently been created in the experiment of Gorlitz et al. [0 . 
This was achieved by lowering the mean-field interaction 
energy in a three-dimensional condensate below the en- 
ergy splitting of either one or two of the directions of 
the harmonic trap, to obtain a two-dimensional or one- 
dimensional condensate, respectively. In a number of 
other experiments a one-dimensional Bose-Einstein con- 
densate was also created in a ^Li-^Li mixture |^ and on 
a microchip 

Theoretically, low-dimensional Bose gases are partic- 
ularly interesting due to the enhanced importance of 
phase fluctuations (^-|8j. Because of these fluctuations, 
Bose-Einstein condensation cannot take place in a ho- 
mogeneous one-dimensional Bose gas at all temperatures 
and in a homogeneous two-dimensional Bose gas at any 
nonzero temperature. This is formalized in the Mermin- 
Wagner-Hohenberg theorem [p|jTc|] . Since this theorem is 
valid only in the thermodynamic limit, it does not ap- 
ply to trapped Bose gases. Therefore, the question arises 
whether under certain conditions we are dealing with a 
true condensate, where the phase is coherent over a dis- 
tance of the order of the size of the system, or only with 
a so-called "quasicondensate" |11 , where the phase is 



coherent over a distance less than the size of the sys- 
tem H-^. This is one of the main questions, which we 



address quantitatively in this paper. 

In the successful Popov theory for three-dimensional 
partially Bose-Einstein condensed gases, the phase fluc- 
tuations are taken into account up to the second order 
around the mean field. In view of the above-mentioned 
importance of phase fluctuations in lower dimensions, 
this is insufficient in general and leads to infrared diver- 
gences. In previous work by three of us, the phase fluc- 
tuations were taken into account exactly [|l2| . The result 
is a mean-field theory which is free of the infrared diver- 
gences in all dimensions. In the present paper, we first re- 
view this modified Popov theory and then extend it to the 
many-body T-matrix theory, by including the effect of 
the medium on the scattering properties of the atoms in 
the gas. The present approach improves on previous at- 
tempts by Petrov et al. to describe low-dimensional 
Bose gases by explicitly incorporating also the effect of 
density fluctuations into the theory. As a result both 
quantum and thermal depletion of the (quasi)condensate 
can now be accounted for and the theory is no longer 
only valid at very low temperatures where the deple- 
tion, and therefore the thermal component in the gas, 
is negligible. Most importantly, we present an equation 
of state for the low-dimensional Bose gas that is free of 
infrared divergences and thus valid in any dimension. For 
a trapped Bose gase this implies that we can determine, 
for a given number of atoms, the density profile of both 
the (quasi)condensate and the thermal cloud in the gas 
for any aspect ratio of the trap. Also the interesting 
crossover problem from a three-dimensional Bose gas to 
a one or two-dimensional one, that is presently being ex- 
plored exprimentally ||l|, can be addressed as well. 

In the present paper we first use the modified many- 
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body T-matrix theory to calculate the one-particle den- 
sity matrix and determine its off-diagonal long-range or- 
der. We also calculate the fractional depletion of the 
(quasi)condensate at zero temperature in one, two, and 
three dimensions. Next, we study the two-dimensional 
homogeneous Bose gas in considerable detail. After hav- 
ing included the phase fluctuations due to vortex pairs 
by a renormalization group approach, we apply the mod- 
ified many-body T-matrix theory to perform an ab ini- 
tio study of the Kosterlitz-Thouless phase transition p^ ] 
from the superfluid to the normal state. Since this is 
a topological phase transition, it cannot be described 
within mean-field theory. Therefore, we proceed as fol- 
lows. We first use the modified many-body T-matrix 
theory to calculate the quasicondensate density and the 
fugacity of vortices. These results are then used as initial 
conditions for a Kosterlitz renormalization group calcula- 
tion. In this manner, we are incorporating critical fluctu- 
ations and able to calculate nonuniversal quantities such 
as the critical temperature for the Kosterlitz-Thouless 
phase transition as a function of the density and the mi- 
croscopic parameters of the theory. 

Finally, we apply the theory to a trapped one- 
dimensional Bose gas, where we calculate the density 
profile at different temperatures. From this we ex- 
tract the crossover temperature for the appearance of 
a (quasi)condensate as a function of the interaction 
strength. We also calculate the behavior of the phase cor- 
relation function which determines whether there exists a 
true condensate or only a quasicondensate. These predic- 
tions are compared to exact results based on a stochastic 
nonlinear field equation for the Bose gas p^ . 

The paper is organized as follows. In Sec. II, we present 
and discuss the Popov theory and its infrared problems. 
We also present our modified mean-field theory in the ho- 
mogeneous limit. In Sec. Ill, we compare the latter with 
exact results in one dimension, and with results obtained 
in the Popov approximation in two and three dimensions. 
We also calculate the reduction of the three-body recom- 
bination rate for a two-dimensional hydrogen gas and 
compare it to the experiment of Safonov et al. |]l5|] . In 
Sec. IV, we study the Kosterlitz-Thouless phase tran- 
sition and in Sec. V, we generalize the many-body T- 
matrix theory to inhomogeneous situations. In partic- 
ular, we consider then a one-dimensional trapped Bose 
gase, for which we compare our predictions to numeri- 
cally exact results. Finally, we conclude and summarize 
in Sec. VI. 



II. MODIFIED POPOV THEORY 

In this section, we derive the modified Popov theory by 
treating the phase fluctuations exactly. We also discuss 
how to incorporate many-body effects into the theory. 
Finally, we give additional arguments for the correctness 



of our approach by using an effective action for the den- 
sity and phase dynamics in a superfluid system that is 
known to give exact results in the long- wavelength limit. 



A. Phase fluctuations 

In order to explain the infrared problems associated 
with the phase fluctuations of the condensate most 
clearly, we first treat a homogeneous Bose gas in a box 
of volume V . Later we generalize to the inhomogeneous 
case. The starting point is the grand-canonical Hamilto- 
nian in second-quantized language 



2™ ^ 



^(x) 



dx / dx>t(x)^t(x') 



xl/(x-x')V'(x')V'(x) , 



(1) 



where /i is the chemical potential, and V^(x) is the atomic 
two-body interaction potential. The mass of the atoms is 
denoted by m, and '(/'^(x) and '(/'(x) are the usual creation 
and annihilation field operators, respectively. 

In the Bose systems considered here and those realized 
in experiment, the temperatures are so low that only s- 
wave scattering is important. Consequently, it is conve- 
nient to neglect the momentum dependence of the inter- 
atomic interaction and use T^(x — x') = Voi5(x — x'). In 
principle this leads to ultraviolet divergences, but these 
can easily be dealt with, as we show later on. The Hamil- 
tonian then reduces to 



2m 



^(x) 



+ -j dxl/o^t(x)^t(x)^(x)^(x) . (2) 

In the presence of a Bose-Einstein condensate the an- 
nihilation operator is parametrized as 



V'(x) 



^'(x). 



(3) 



where uq is the condensate density and '0'(x) describes 
the fluctuations. The standard one-loop expressions for 
the density n and the chemical potential /i are obtained 
after a quadratic approximation to the Hamiltonian in 
Eq. (^, i.e., by neglecting terms that are of third and 
fourth order in the fluctuations. This yields [0,^61 



Ck + mVo 



(4) 
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, 2ek + noVb , ^ 

TlLOVr 



(5) 



where ?iwk = (^k + 2fioK)ek)^^^ is the Bogohubov disper- 
sion relation, N{x) — XjieP'-^ — 1) is the Bose-Einstein 
distribution function, and /? — 1/kBT is the inverse ther- 
mal energy. 

In agreement with the Mermin-Wagner-Hohenberg 
theorem, the momentum sums in Eqs. (|) and (|) con- 
tain terms that are infrared divergent at all temperatures 
in one dimension and at any nonzero temperature in two 
dimensions. The physical reason for these "dangerous" 
terms is that the above expressions have been derived by 
taking into account only quadratic fluctuations around 
the classical result no, i.e., by writing the annihilation 
operator for the atoms as ■(/'(x) = + V''(x) and ne- 
glecting in the Hamiltonian terms of third and fourth 
order in As a result the phase fluctuations of the 

condensate give the quadratic contribution no(x(x)x(x)) 
to the right-hand side of the above equations, whereas an 
exact approach that sums up all the higher-order terms 
in the expansion would clearly give no contribution at 
all to these local quantities because no(e~'*(''^e**^''^) = 
"0 (1 + (x(x)x(x)) -I- . . .) = 1. To correct for this we thus 
need to subtract the quadratic contribution of the phase 
fluctuations, which from Eqs. ^ and (||) is seen to be 
given by 



"o(xWx(x)) 



-E 

V ^ 



2?ia;k 



[1 + 2N{hoJi,)] 



(6) 



As expected, the infrared divergences that occur in the 
one and two-dimensional cases are removed by perform- 
ing this subtraction. 

After having removed the spurious contributions from 
the phase fluctuations of the condensate, the resulting 
expressions turn out to be ultraviolet divergent. These 
divergences are removed by the standard renormalization 
of the bare coupling constant Vq ■ Apart from a subtrac- 
tion, this essentially amounts to replacing everywhere the 
bare two-body potential Vb by the two-body T-matrix 
evaluated at zero initial and final relative momenta and 
at the energy —2//, which we denote from now on by 
T2B(--2^). It is formally defined by 
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r2B(-2/i) 



1 1 



1 



1/ ^ 2ek 

k 



2^ 



(7) 



Note that the energy argument of the T-matrix is —2/^, 
because this is precisely the energy it costs to excite two 
atoms from the condensate ||l^,^. After renormaliza- 
tion, the density and chemical potential are 



no 



-E 
V ^ 



ek 



2huo\, 



nor2B(„2M) 
2ek + 2/i 



ek 

flUJv, 



11= (2n-no)r2"(-2^) 
= (2n' + no)r2B(-2^) , 



(8) 



(9) 



where n' = n — no represents the depletion of the conden- 
sate due to quantum and thermal fluctuations and the 
Bogoliubov quasiparticle dispersion now equals huj]^ = 

1/2 

[e^ -I- 2nor2^(— 2/x)ek] . The most important feature 
of Eqs. (||) and (^) is that they contain no infrared and ul- 
traviolet divergences and therefore can be applied in any 
dimension and at all temperatures, even if no condensate 
exists. 

Note that Eq. (||) is also ultraviolet divergent. The 
ultraviolet divergences are removed by the same renor- 
malization of the bare interaction Vq and the final result 



(x(x)x(x)) = 



r-^"(-2/x) 

V 

1 

2ek + 2fj, 



E 



2hiLj]i 



[1 + 2N{hujk)] 



(10) 



We will return to the physics of this important expression 



in Sec. II C below 



B. Many-body T-matrix 

In the previous section, we presented the modified 
Popov theory that takes the phase fluctuations into ac- 
count exactly. The final results in Eqs. (||), and ( p^ ) 
involve the two-body T-matrix. The two-body T-matrix 
takes into account successive two-body scattering pro- 
cesses in vacuum. However, it neglects the many-body 
effects of the surrounding gas. In order to take this into 
account as well, we must use the many-body T-matrix in- 
stead of the two-body T-matrix in Eqs. (||), @ and 
Many-body effects have been shown to be appreciable in 
three dimensions only very close to the transition tem- 
perature [l9|] , but turn out to be much more important 
in one and two dimensions Since the effects of the 
medium on the scattering properties of the atoms is only 
important at relatively high temperatures, we can apply 
a Hartree-Fock approximation to obtain for the many- 
body T-matrix 



r'^^"(-2/x) =T^^(-2Ai) 



1 + T2B(_2^)^^ 



Niek + noT^^i-2^i)) 

Ek + M 



(11) 
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The situation is in fact slightly more complicated be- 
cause we actually have two coupling constants in the 
equation for the chemical potential, which is the homo- 
geneous version of the Gross-Pitaevskii equation. When 
two atoms in the condensate collide at zero momentum, 
they both require an energy fi to be excited from the con- 
densate, and thus the coupling is evaluated at — 2/x. This 
is the coupling that multiplies uq in the Gross-Pitaevskii 
equation. On the other hand, the coupling that mul- 
tiplies n' in the Gross-Pitaevskii equation involves one 
condensate atom and one atom in the thermal cloud, so 
that this coupling should now be evaluated at — /i. The 
equation for the chemical potential thus becomes 



(12) 



Note that the existence of two different many-body cou- 
pling constants for the interatomic interactions has pre- 
viously been discussed by Proukakis et al. . (See how- 
ever also Ref. This lead these authors to the so- 
called Gl-theory, which is qualitatively somewhat similar 
to Eq. (p^ ) but differs in detail. 



C. Long- wavelength physics 

We have given physical arguments for how to identify 
and subtract the contribution to Eqs. (Q) and (^) from 
the phase fluctuations of the condensate. At this point, 
we would like to give a somewhat more rigorous field- 
theoretical argument. The Euclidean action that corre- 
sponds to the Hamiltonian in Eq. (0) is 





1 rt^P 



d 



dr 







dr 



2m 



4, 



(13) 



If we substitute -(/"(x, t) = \Jn + Sn{x, T)e^^'-^''^^ into 
Eq. (|l3|), we obtain the action 



S[Sn, x] = 



hp 



dr / c?x 



1 ^ dSn 
2^57 



-ih{n + Sn)^ + {Vxf 
f 

-—on 

2 \ 2mn 



+ Sn 



(14) 



Here, n is the average total density of the gas and 
(5ri(x, r) represents the fluctuations. At zero tempera- 
ture, this action is exact in the long-wavelength limit, if 
(fi^k'^ /2mn + Vq) is replaced by Xm(k), where Xnn(k) is 
the exact static density-density correlation function. 

By using the classical equation of motion to eliminate 
the phase x(^? obtain the following action for the 

density fluctuations 5n(x, r): 



i-np r 
S[6n] = y dr J 



. dSr. 



dr 



, dSn 



+ 2^'^X„n(-*V)(5n 



(15) 



The density fluctuations are therefore determined by 

1 



(a(x)a(x')) = ^ E ^ 



xe 



!k-(> 



(16) 

where w„ = 27rn/h(3 are the even Matsubara frequencies 
and TiWk = Vn^k/XrmOO • Summing over these Matsub- 
ara frequencies, we obtain 

(5n(x)5n(x')) = ^ E ^ [1 + ^Nihv^)] 

k 



^ik-(x-x') 



(17) 



Similarly, by using the classical equation of motion for 
(5ri(x, r), we obtain from Eq. (|lj) the following action for 
the phase fluctuations: 



Six] 



hp 



dr i dx 



.(-iV 



2m ^ -^^ 



(18) 



From this action, it is straightforward to calculate the 
propagator for the field xk^-, ''") ^^id thereby the correla- 
tion function (x(x)x(x')). The result is 

^ ^ [1 + 2N{riujy,)\ 



(X(x)x(x')) = -^E 



^ Xnn 
■ik- (x— x') 



(k) 2huj]^ 



xe 



(19) 

Setting x' = x, we recover Eq. (nffl) in the long- 
wavelength limit, if we use Xnn(k) ~ 1/T^^(— 2/i) for the 
static density-density correlation function in that limit. 
It is important to mention that Eq. (^9|) is often used for 
the short-wavelength part of the phase fluctuations as 
well [0,1). This is, however, incorrect because it contains 
ultraviolet divergences due to the fact that the above pro- 
cedure neglects interaction terms between density and 
phase fluctuations that are only irrelevant at large wave- 
lengths. The correct short wavelength behavior is given 
in Eq. m. 



III. COMPARISON WITH POPOV THEORY 

We proceed to compare predictions based on Eqs. (|^), 
(^, and ( p^ ) with exact results in one dimension and 
results based on the Popov theory in two and three di- 
mensions. We consider only the homogeneous case here 
and discuss the inhomogeneous Bose gas in Sec. 
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A. One dimension 

To understand the physical meaning of the quantity no 
in Eqs. (||) and (joj), i.e., whether it is the quasiconden- 
sate density or the true condensate density, we must de- 
termine the ofF-diagonal long-range behavior of the one- 
particle density matrix. Because this is a nonlocal prop- 
erty of the Bose gas, the phase fluctuations contribute 
and we find in the large- |x| limit 



(20) 



Using Eq. ([To|), we obtain for the exponent in Eq. (pO|) 



([x(x)-x(or) 



V 



E 



1 



[1 + 2N{huJk)] 



k L 

X [1 — cos(k-x)] 



(21) 



Writing the sum over wave vectors k as an integral, the 
phase fluctuations at zero temperature can be written as 



iim-mn 



1 



27rno^ 



1 — cos(fca;) 



dk 

k\/k'^ 



1 



(22) 



where $, = h/[AmnoT^^{—2^)Y/^ is the correlation 
length. Note that we have used that r^^(-2^) = 
T^^(— 2/i) at zero temperature and that the chemical po- 
tential, as we show shortly, is to a good approximation 
equal to noT^^(— 2/z). The integration can be performed 
analytically and the result is 



[X(x)-x(0)]') = 



1 



27rnoC 



^1^2(1/2; 1, 3/2; 



2^2 



2^3(1,1; 3/2, 3/2; 2x74^2) 



(23) 



where iFj{ai, a2, -..ai; Pi, P2, ■■■Pj; x) are hypergeometric 
functions. In the limit |x| cx), Eq. (|23) reduces to 



([X(x)-X(0)]') 



1 



27moC 



log(x/^) 



(24) 



Using Eq. (^J) we find that the one-particle density ma- 
trix behaves for Ixl — )■ cx3 as 



(V't(x)^(O)) 



HQ 



(x/e)V 



(25) 



A few remarks are in order. First, the asymptotic be- 
havior of the one-particle density matrix at zero temper- 
ature proves that the gas is not Bose-Einstein condensed 
and that no should be identified with the quasiconden- 
sate density. Second, in the weakly-interacting limit 



inn^ ^ 1 the depletion is small, so that, to first approxi- 
mation, we can use ng ~ n in the exponent 77 = l/AnnQ^. 
Indeed, from Eqs. (||) and ^ we obtain the following 
expression for the fractional depletion of the quasicon- 
densate 



n — riQ 



n 



■^/2 



47rn^ 



TT - 1 



(26) 



We see that the expansion parameter is l/47rn^ and, 
therefore, that the depletion is very small. Keeping this 
in mind, Eq. (|2^) is in complete agreement with the ex- 
act result obtained by Haldane . Note that our theory 
cannot describe the strongly- interacting case Aim^ <C 1, 
where the one-dimensional Bose gas behaves as a Tonks 
gas 

Finally, our results show that at a nonzero temperature 
the phase fiuctuations increase as ([x(x) — x(0)]^) cx: |x| 
for large distances, and thus that the off-diagonal one- 
particle density matrix vanishes exponentially. Hence, at 
nonzero temperatures not even a quasicondensate exists 
and we have to use the equation of state for the normal 
state to describe the gas, i.e., 



(27) 



where the Hartree-Fock self-energy satisfies 

2nT^^(-?iE ) , (28) 
and the many-body T-matrix obeys 

T^^i-hT.) = r2B(_?iE) 

iV(ek + TiE - Ai)) 



Ek + ?iS/2 



(29) 



Note that the last three equations for the description of 
the normal phase of the Bose gas are again valid for an 
arbitrary number of dimensions. 



B. Two dimensions 



In analogy with Eq. (^2|) , we obtain for the phase fiuc- 
tuations in two dimensions at zero temperature 



[X(x)-X(0)]'> 



1 



^dkl^^, (30) 



2n^noC Jo Vk^ + 1 

Therefore we now find in the limit |x| — )■ 00 that 

(7^t(x)v3(0)) =no 



(31) 



and no is clearly the condensate density of the gas. How- 
ever, at nonzero temperatures the correlation function 
behaves as 
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(^t(x)^(O)) 



no 



l/noA= 



(32) 



where A = y / rak-oT is the thermal de Broghe wave- 
length and no corresponds again to the quasicondensate 
density. At zero temperature, the fractional depletion 
of the condensate in the Popov approximation was first 
calculated by Schick IT^]. He obtained 



^ = f T-(-2,) 



(33) 



where the chemical potential satisfies fi = nT^^(— 2/i). 
The corresponding result based on Eqs. (||) and (O) is 



^^-f (l-ln2)T2^(-2^) 
n 47r 



(34) 



where /i now satisfies Eq. (|^). In two dimensions, the de- 
pletion predicted by the Popov theory is thus too large 
by a factor of approximately 3. 

In a number of applications, we need to calcu- 
late many-body correlators. For instance, in or- 
der to calculate how a quasicondensate modifies the 
two-body relaxation constants of a spin-polarized two- 
dimensional Bose gas, we need to know K<-'^^{T) = 
(^t(x)?/''l'(x)?/;(x)?/;(x))/2n^. This correlator was consid- 
ered in Ref. ||25|| using the many-body T-matrix theory 
with an appropriate cutoff to remove the infrared diver- 
gences. An exact treatment of the phase fluctuations 
leads however directly to an infrared finite result as we 
show now. Using the parametrization in Eq. (^) for the 
annihilation operators, we obtain first of all 

(^t(x)^t(x)^(x)^(x)) =nl + no [(^'(x)^'(x)) 

+ (^'t(x)^'t(x)) + 4(V>'t(x)V;'(x))" 
+2(^'t(x)^'(x))2 

+ (^'(x)^'(x))(V3't(x)Vi't(x)) . (35) 

The normal average is given by (■!/''^(x)'0'(x)) = 
'T-' + ^o(x(x)x(x)) and the anomalous average obeys 
(■0'(x)^/;'(x]) = — rio(x(x)x(x)), as we have seen. Using 
this, Eq. ( |35| ) can then be written as 

(^t(x)V?(x)Vi(x)Vi(x)) ^nl[l + 2(x(x)x(x)) 
+3(x(x)x(x))2] 

+4no [1 + (x(x)x(x))] n + 2{n'f . (36) 

Writing the correlator in this form, we explicitly see that 
the infrared divergences are due to spurious contributions 
from the phase fluctuations. Removing them, we obtain 
for the renormalized correlator 



1 
2^ 



'\2 



^o + 4r7.on' + 2 (n') 



(37) 



We would like to point out that critical fluctuations 
are not treated within our mean-field theory. This is of 
course essential in the study of the Kosterlitz-Thoulcss 



phase transition and we return to this issue in Sec. IV. 
However, an example of a physical observable where 
phase fluctuations are not important, is the three-body 
recombination rate constant. We are at this point, there- 
fore, already in the position to determine the reduction 
of the three-body recombination rate constant due to the 
presence of a quasicondensate. This can be expressed 



L{T) 



r2B(-2?iE) 



(38) 



where is the recombination rate constant in the nor- 
mal phase, which is essentially independent of temper- 
ature, and the self-energy satisfies = 2nT'^^ {-hJ:) . 
The renormalized three-body correlator 



9nln' 



187io(7i')^ + 6(n')^ (39) 



is obtained from the expression for the correlation func- 
tion (V'^(x)?/;^^(x)'!/;t(x)?/)(x)?/;(x)?/;(x)) by removing, as 
before, the spurious contributions from the phase fluc- 
tuations. Moreover, in two dimensions the T-matrix de- 
pends logarithmically on the chemical potential as 



T^^{-2pi) = 



47rn^ 



I 



m h\{2n^ / ^ima^) 



(40) 



where a is the two-dimensional s-wave scattering length. 
In the case of atomic hydrogen adsorbed on a super- 
fluid helium fllm, the scattering length was found to be 
a — 2.4ao p6|| , where oo is the Bohr radius. However, 
there is some uncertainty in this number because the hy- 
drogen wave function perpendicular to the helium surface 
is not known very accurately. In order to compare with 
experiment, we may therefore allow a to vary somewhat. 

In Fig. ^ we show the reduction of the three-body re- 
combination rate as a function of the density at a fixed 
temperature T = 190 mK for three different values of 
a. As can clearly be seen from Fig. ^ the reduction of 
the three-body recombination rate is very sensitive to the 
value of a. What is most important at this point is that at 
high densities our calculation shows that the reduction of 
the recombination rate is much larger than the factor of 6 
predicted by Kagan et al. . Such large reduction rates 
are indeed observed experimentally Jist . A direct com- 
parison, however, between the results of our theory and 
the measurements of Safonov et al. cannot be made here, 
since the density and temperature of the adsorbed hydro- 
gen gas were not measured directly, but inferred from the 
properties of the three-dimensional buffer gas. Because 
this procedure requires knowledge of the equation of state 
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of the two-dimensional Bose gas absorbed on the super- 
fluid hehum film, the raw experimental data needs to be 
reanalyzed with the theory presented in this paper. We 
can, however, compare the density at which the recombi- 
nation rate starts to deviate considerably from the result 
in the normal state. For the temperature of T = 190 
mK, where most of the experimental data is taken, this 
is at a density of about 1.0 x 10^'^ cm"^, which is in excel- 
lent agreement with experiment. In view of this and the 
above mentioned problems we thus conclude that our re- 
sults present a compelling theoretical explanation of the 
experimental findings. 



C. Three dimensions 

The Popov theory has been very successful in describ- 
ing the properties of dilute three-dimensional trapped 
Bose gases. It is therefore important to check that an ex- 
act treatment of the phase fluctuations leads at most to 
small changes in the predictions for the three-dimensional 
case. 

At zero temperature, the fractional depletion within 
the Popov theory was first calculated by Lee and 
Yang pi] and is given by 



n — no o / 7ia 



n 6 \ TT 

where a is the s-wave scattering length and we have used 



Anah 



The result that follows from Eqs. (^) and (||) is 



(42) 



(43) 



The fractional depletion is approximately 2 /3 of the value 
obtained from the Popov theory. It turns out that this 
is the largest change in the condensate depletion, since 
the effects of phase fiuctuations decrease with increasing 
temperature. The critical temperature Tbec is found by 
taking the limit no ^ in Eqs. (||) and (^. These expres- 
sions then reduce to the same expressions for the density 
and chemical potential as in the Popov theory. This im- 
plies that our critical temperature for Bose-Einstein con- 
densation coincides with that obtained in Popov theory, 
i.e., the ideal gas result 



Tbec — 



27rr 
mkn 



c(f] 



2/3 



(44) 



IV. KOSTERLITZ-THOULESS PHASE 
TRANSITION 

In the previous section, we have compared our results 
using the modified many-body T-matrix theory with es- 
tablished results in one, two, and three dimensions in the 
Popov approximation. Due to the mean-field nature of 
the modified many-body T-matrix theory, the Kosterlitz- 
Thouless transition is absent and a nontrivial solution of 
the equation of state exists even if the superfiuid den- 
sity Us obeys risA^ < 4. In this section, we correct for 
this by explicitly including the effects of vortex pairs in 
the phase fluctuations. The idea is to use the modified 
many-body T-matrix theory to determine the initial val- 
ues of the superfiuid density and the vortex fugacity, and 
to carry out a renormalization-group calculation to find 
the fully renormalized values of these quantities. In this 
manner we can for example calculate the critical temper- 
ature Tc for the Kosterlitz-Thouless transition given the 
scattering length a and density n. 

Let us for completeness first briefly sketch the deriva- 
tion of the renormalization group equations for the su- 
perfluid density and the vortex fugacity |2|]. Consider 
the velocity field of a vortex where the core is centered 
at the positions x^, which we for simplicity take to lie 
on a lattice with an area of the unit cell equal to f2. By 
rotating the velocity field by ninety degrees we can map 
it onto the electric field of a point charge in two dimen- 
sions. Since the total energy in both systems is propor- 
tional to the square of the field integrated over space, 
there is complete analogy between a system of vortices 
and a two-dimensional Coulomb gas. This analogy is very 
useful and we will take advantage of it in the following. 
The total vorticity corresponds to the total charge of the 
Coulomb gas. For the analogous two-dimensional neutral 
Coulomb gas on a square lattice, the partition function 
can be written as 



(45) 



where V^(xi — Xj) — — 27r?i^ns ln(|xi — Xj|/^)/m is the 
Coulomb interaction between two unit point charges in 
two dimensions, is the superfiuid density, and Ec is 
the energy associated with the spontaneous creation of a 
charge, i.e., it is the core energy of the vortices. The sum- 
mation is over all possible conflgurations of charges Ui at 
positions x.^ on the lattice. The partition function can 
be rewritten in a field-theoretic fashion in terms of the 
electrostatic potential 0(x) and the fugacity y = e~^^'' 
as 



where C, 



2.612. 



z= y ^<^^ 



(46) 
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where K' = In the Hmit where y <C 1, 

the charge density is very low, and thus only rij = 0, ±1 
contribute to the partition function. We can then write 

I l?0e-/''"i^'(^*)'[] [l + yexp(z/30(xj)) 
•' j 

+yexp(-i/30(xj)) + ...j 

Vcf) e~/ d^[^K'{V4>?-gcos(l34,)\ ^ ^^^^ 



where g — 2y /fl. It is convenient to introduce a dimcn- 
sionless dielectric constant K that is related to K' by 
K = P'^/Att'^K' = nsA'^/2'K where A is the thermal wave- 
length. 

The renormalization group equations for K, which 
is thus proportional to the superfluid density, and the 
fugacity y can now be obtained by performing the 
usual momentum-shell integrations. For the Sine-Gordon 
model derived in Eq. (M^), this results in 



dK-^{l) 
dl 

dyjl) 

dl 



= 47rV(0 + O(2/3) , 

= [2-nK{l)]y{l) + 0iy') . 



(48) 
(49) 



The renormalization group equations to leading order in 
the variables K{1) and y{l) were first obtained by Koster- 
litz |3^], while the next-to- leading order terms were de- 
rived by Amit et al. The flow equations are not 
significantly changed by including the higher-order cor- 
rections and we do not include them in the following. 

The renormalization group equations (^ ) and ( ^ ) can 
be solved analytically by separation of variables and the 
solution is 



I 

2^ 



K{1) 



n\og{K{l)) 



(50) 



where the integration constant C is determined by the 
initial conditions. For the critical trajectory it can be 
calculated by evaluating the left-hand side at the fixed 
point {y{oo), K(oo)) = (0, 2/7r). In this manner, we find 
C = [log(7r/2) - l]/27r2 ~ --0.0278. In Fig |, we show the 
flow of the Kosterlitz renormalization group equations. 
There is a line of fixed points y(oo) = and K{oo) > 0. 
The fixed point {y{(X)), K{oo)) — {0,2/tt) corresponds to 
the critical condition for the Kosterlitz-Thouless transi- 
tion, where the vortices start to unbind and superfluidity 
disappears. Physically this can be understood from the 
fact that below the transition the fugacity renormalizes 
to zero, which implies that at the largest length scales 
single vortices cannot be created by thermal fluctuations. 
They are therefore forced to occur in pairs. 

The initial conditions for the renormalization group 
equations are 



K{0) 

y(0) 



t;2 

n no 
mk^T 



(51) 
(52) 



where no is the quasicondensate density and Ec is the 
core energy of a vortex. Both are obtained from the mod- 
ified many-body T-matrix theory considered previously. 
Writing the order parameter for a vortex configuration 
as "00 (x) = v^^/(a;/^)e*'', where i? is the aziniuthal an- 
gle, the core energy of a vortex follows from the Gross- 
Pitaevskii energy functional. It reads 



Ec = -z—no'n: I dx 
Zra 



2lf 

. dx 



(53) 



The dimensionless integral was evaluated by Minnhagen 
and Nylen, and takes the value 1.56 

Using the solution to the flow equations ( |50| ) and the 
initial conditions, we can calculate the temperature for 
the Kosterlitz-Thouless transition given the scattering 
length a and the density of the system. In the follow- 
ing, we consider again atomic hydrogen. In Fig. ^ we 
show the critical temperature as a function of density for 
a = 2.4ao. We see that the critical temperature is essen- 
tially proportional to the density of the system. This can 
be seen in more detail in Fig. ^, where we plot nA^ as a 
function of n. It is clear from this figure that nA^ indeed 
changes only slightly over the density range considered. 

To understand the physics of the calculation better, we 
show in Fig. |^ the quasicondensate fraction no/n follow- 
ing from the many-body T-matrix theory as a function of 
temperature for a total density n — 1.25 x fO^'^cm^^. In 
addition, we show the superfluid density tIj as calculated 
from the renormalization-group procedure explained pre- 
viously. The Kosterlitz-Thouless transition takes place 
when ris lies on the line given by ngh? = 4. Noticing 
that the left-hand side of Eq. (^0|) is a function of uqA^ 
only and solving the equation with respect to noA^ using 
the value of C at the transition, we obtain the condition 
ngA^ ~ 6.65 for the Kosterlitz-Thouless transition. It is 
therefore also seen in Fig. || that the Kosterlitz-Thouless 
transition takes place when the line given by tiqA^ ~ 6.65 
intersects with the curve for ng- 



V. TRAPPED BOSE GASES 



In this section, we generalize the theory presented in 
Sees. ^ and ID to inhomogeneous Bose gases. We also 
apply the results to a trapped one-dimensional Bose gas. 
We start by generalizing our previous expressions for 
the total density, Eq. (H), and the phase fluctuations, 
Eq. (p^, to the inhomogeneous case. To do so we first 
consider the Gross-Pitaevskii equation 
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— + y^^'(x) + T^^(-2a*(x))|V'o(x)P 



+ 2T^^(-A*(x))n'(x) 



■00 (x) = mV'o(x) 



(54) 



which generahzes Eq. ^ to trapped Bose condensates. 
Here the local chemical potential equals /i(x) = fj, — 
V°^*{x.). The noncondensed density n'(x) is to be de- 
termined by solving the Bogoliubov-de Gennes equations 



"i(x) 



+ T^^^(-2M(x))no(x)«,(x) , 



-huj-iVj{:>c) 



Wj(x) 



(-2^(x))no(x)Mj(x) , 



(55) 



(56) 



where no(x) = |'(/jo(x)p and the Hartree-Fock potential 
V^^{x.) is given by 



^"^(x) = V^^'^'(x) + 2T^^(~A*(x))n'(x) 



(57) 



The functions Uj and Vj are the usual Bogoliubov particle 
and hole amplitudes respectively, which are chosen to be 
real here. In some cases, for instance when ■00 describes 
a vortex, we cannot choose these amplitudes real and our 
equations are easily generalized to incorporate that fact. 

In terms of the Bogoliubov amplitudes, the expression 
for the total density in Eq. (S) reads 



i(x) = no(x) +^ 



+ ^'j(x)(wj(x) +Uj(x)) 

^ TMB(_2^(x))no(x) 
2ej + 2/i(x) 



(0.(x))' 



(58) 



Here, (f>j is the large-j or high-energy limit of Uj which 
can be obtained by neglecting the interaction terms in 
Eq. (p5|), namely 



(59) 



= ( -^v^ + y-*(x) - M ) 0,(x) 



In the large-j limit, we also have 

TMB(-2M(x))no(x) 



2ei 



(/)j(x) 



(60) 



It is clear that the expression of Eq. ( |5q ) for the total 
density is ultraviolet finite since the second and third 
term cancel each other in the large-j limit. 

Finally, the phase fluctuations in the trapped case are 
determined by (x(x)x(x')) which is given by 



(X(x)x(x')) 

-E 



1 



J 2\/no(x)no(x') \ 

r^^(_2^(x))no(x) 
2ej + 2/i(x) 



Uj{x')vj{x) [l + 2N{hujj)] 



0,(x')0,(x) 



+Uj{x)vj{x') [I + 2N{hLjj)] 



r^"(-2A.(xO)no(xO 
2ej + 2^(x') 



0j(x)(/'j(x') 



(61) 



In particular the normalized form of the off-diagonal one- 
particle density matrix of Eq. (20) becomes for large dis- 
tances |x — x'l equal to 



5(i)(x,x') = exp(-([x(x)-x(x')n/2) 



(62) 



A. Density profiles 

We are now ready to calculate the total density pro- 
file by solving Eqs. ( |5^ ) and (|5^) selfconsistently. In 
the remainder of the paper, we restrict ourselves to one- 
dimensional harmonic traps with, therefore. 



1 



2 -,2 



(z) = —VTLLOZZ 
^ ' 2 



(63) 



For simplicity we use the local-density approximation, 
which allows us to calculate the densities directly using 
the many-body generalization of Eq. (^, and Eq. (p^). 
In Fig. H the total density profile is shown at four dif- 
ferent values of the temperature. For the four different 
temperatures each of the four curves is composed of two 
parts. The first part near the center of the trap rep- 
resents the superfluid part of the gas and contains the 
(quasi)condensate. The other part consists only of the 
noncondensed atoms. The small discontinuity between 
the two parts is caused by the use of two different equa- 
tions of state for the superfluid and thermal phases of the 
gas. In the following we call the position of the disconti- 
nuity the temperature-dependent Thomas-Fermi radius 
of the (quasi)condensate. For distances below the dis- 
continuity we use the above-mentioned equations, while 
for distances above the discontinuity, we simply use 



f°° dk 

7j{z)= / — N{ek + mLulz^/2 
J-oo 27r 

^^'^i~hJ:iz))-^,) . 



2nT" 



(64) 



For all these curves fi = SOfkUz- The remaining parame- 
ters used here are those of the experiment of Gorlitz et 
al. [Q. In particular, we have used ^^Na in the trap with 
ujz = 2iT X 3.5 rad/sec, Iz = ^JfijmbJz — 1.12 x 10~^m. 
The three-dimensional s-wave scattering length is a ~ 
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2.75 nm which is related to the one-dimensional scatter- 
ing length defined by r^^(— 2^) = A'KKh'^/m. For 
harmonic confinement, we have k = a/2'Kl\, where 
is the harmonic oscillator length of the axially symmet- 
ric trap in the direction perpendicular to the z-axis. We 
have used uj±^ = 27r x 360 rad/sec and = ^Jh/nuo±^ ~ 
1.10 X lO-'^m. 

As expected, the temperature-dependent Thomas- 
Fermi radius decreases with increasing the temperature. 
At the temperature when this radius vanishes the one- 
dimensional system reaches the crossover temperature for 
the formation of a (quasi)condensate. We have calculated 
this crossover temperature for different values of the scat- 
tering length at a constant value of the number of atoms, 
the latter being fixed by adjusting the chemical poten- 
tial. In Fig. 1^, we show the result of this calculation, and 
plot the crossover temperature Tqc and the chemical po- 
tential against the scattering length. The inset in Fig. 
shows that on a double logarithmic scale the tempera- 
ture Tqc is clearly not a straight line indicating that the 
relation between Tqc and k is not a simple power law 
and may contain logarithmic dependence. It is shown in 
Ref. that for a — the transition temperature Tqc 
should satisfy Tqc = NhLLjz/kB^n{2N), where N is the 
number of atoms. In the case of Fig. |^ we have N = 950, 
which leads to Tqc — 164To for an ideal gas. Of course, 
this limit is not obtained in Fg. |^ because our calcu- 
lation is based on a local-density approximation, which 
will always break down for sufficiently small values of n. 
On the other hand, the curve for the chemical potential 
becomes almost a straight line on a double logarithmic 
scale. A calculation of the slope of this line shows that 
the slope starts at a value slightly larger than 2/3 at 
the lower end of the curve and saturates at this value 
near the upper end. The value 2/3 is what we expect, 
since in the Thomas-Fermi limit it is easy to show that 
H = {3Tr/V2f/^{NK)'^/^huj, ~ 3.5{NK)^/^nuj,. Calculat- 
ing similar curves for different values of N we actually 
find numerically that fj, ~ 3.2{Nn)'^/'^huJz- 



B. Phase fluctuations 

The aim of this section is to calculate the normal- 
ized off-diagonal density matrix given by Eq. (|6^). This 
function expresses the coherence in the system. It is 
calculated by solving the Bogoliubov-de Gennes equa- 
tions in Eqs. (^5|) and (|5^) using the density profile cal- 
culated in the previous subsection. Specifically from 
the (quasi)condensate density profile uq, we determine a 
temperature-dependent Thomas- Fermi radius. This ra- 
dius is then used to calculate the phase fluctuations at 
that specified temperature in the following manner. 

We start by employing the following scaling: lengths 
are scaled to the trap length Iz = {Ti/mujz)^^'^ , frequen- 
cies to Wz, energies to TitOz, and densities to 4:Tt/Iz- With 



this scaling, the Bogoliubov-de Gennes equations take 
the dimensionless form 



1 



1 



Lo-iUi = I ^2^^ ^2^ ^ M + 2Kn ) Uj — kuoVj , (65) 



Id^ 
2d^ 



+ -Z — IJ.+ 2Kn I Vj — KTlQUj 



(66) 



Using the same scaling, the Gross-Pitaevskii equation 
takes the form 



1 d- 



1 



no = 



(67) 



Next we define Fj{z) = Ujlz) + Vj{z) and GAz) ~ 
Uj{z) — Vj{z), and derive from Eqs. (65) and ( pq ) two 
equations for Fj{z) and Gj{z), namely 



d^F 
1^ 



2 ^dgdF 
dz^ dz dz 



d^G d^G df dG 

-2(/ + 5)^-4^^ 



2g-4./)G = 



where the functions f{z) and g{z) are given by 
1 



2 

1 , 



2Kn — n 
2Kn 



KUq , 



/i — KUq . 



(68) 

(69) 

(70) 
(71) 



For our purposes we can use the Thomas-Fermi approx- 
imation which neglects the derivative term in Eq. ( |67 
Hence 



— /i + K{no + 2n') 



= . 



(72) 



In this limit, the functions f{z) and g{z) are given by 
f{z) — 2Kno(z) and g{z) = 0. In the Thomas-Fermi 
approximation, we substitute these values for f(z) and 
g{z) into Eqs. ( |6^ ) and (|69| ) and neglect the fourth-order 
derivative terms. These Equations thus take the form 



dFf 



0, 



d'^lnnoGj) 
d^ 



.|G, = . 



(73) 
(74) 



In Ref. |3J|, it was shown that y/ KnQ{z)Gj{z) corre- 
sponds to density fluctuations and Fj{z) / nriQiz) cor- 
responds to phase fluctuations in the hydrodynamic ap- 
proach We therefore define the function hj{z) 
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(75) 



where 



Substituting this back in Eqs. ( |73| ) and ( |74| ) both equa- 
tions reduce to a single equation for hj{z), namely 



(Ph^ dnn dhi ^, 

^ ' K-r--r^ +UJih., = 

dz dz ■> 



dz^ 



(76) 



This equation can finally be simplified using the Thomas- 
Fermi expression for Kno{z) from Eq. ([72[), namely 
Kno(z) ~ /i' — where ^' = ^ — 2Kn'(0). Note that 
we have made the approximation that we take n'{z) to 
be equal to its value at the center, namely n'(0). This 
approximation is justified in view of the fact that the 
presence of the condensate repels the atoms from the 
noncondensate atoms from the center of the trap. This 
is also supported by a numerical solution of Eqs. (||) 
and (p^), where we find that n'(z) <C ^0(2;), except at 
the Thomas- Fermi radius where they become of the same 
order. Moreover, the slope of n'{z) is small for distances 
close to the center. Thus, the last equation becomes 

(1 - y')^^^(y) - 2y^/»,(y) + 2^]h,iy) = , (77) 



where y = z/Rtp{T) and Rtf{T) = is the 

Thomas-Fermi radius 

In the following, we reinstate the units. Interestingly, 
Eq. ( [77] ) is the Legendre equation with the Legendre poly- 
nomials as solutions: 

h,{z)^P,{z/RTF)^Pj{y), (78) 
where the energy eigenvalues are 



tiuji 



jU + 1) 



J = 0,l,2,. 



(79) 



The normalization condition for the Bogoliubov ampli- 
tudes is 



R.T 

-Rtf 

which leads to 



dz (z)-|«,f (z)] =1, 



(80) 



\ARtf V ri^^] 



Gj{z) 



1 {3 + 1/2)?^ J P,{y) 



TF 



(y) , (81) 



(82) 



These expressions are in agreement with those obtained 
in Ref. Bfl]. Consequently, we find 




1 



'(j + 1/2)^' 



\ARtf y 

1 Uj + l/2)hiuj 



TF 



(85) 
(86) 



The expression for the phase fluctuations in Eq. ( |6l| ) now 
reads, after neglect of the quantum contribution. 



([X(z)-x(^')] > = 



n2 

"tf 



J2N{nu;,)\A'Jp,{y)~P,{y' 



1 2 



I — y^ I — y'- 



(87) 



It should be noted that the first term in this sum, j = 0, 
does not diverge as one might think at first instance. 
It actually vanishes and the sum can start from j = 1. 
Physically, this is a result of the fact that the global phase 
does not influence the phase fluctuations. 

For the four values of temperature used in Fig. ^, we 
insert the corresponding Rtf{T) in Eq. ( |87| ) to calculate 
the phase correlation function (;^^^(0, z). In Fig. ||, we 
plot this quantity and we see that at sufficiently low tem- 
peratures the phase correlation function decreases only 
slightly over the condensate size. This indicates that a 
true condensate can exist at sufficiently low temperatures 
for interacting trapped one-dimensional Bose gases. 



C. Comparison with exact results 



We next compare the above results to predictions 
based on a Langevin field equation for the order param- 
eter of a trapped, one-dimensional condensate in con- 
tact with a three-dimensional Bose gase that acts as a 
"heat bath" . Such a situation can be created experimen- 
tally in a magnetically trapped three-dimensional sys- 
tem, by using a laser beam to provide and additional 
optical potential along two of the directions. The laser 
beam then needs to be focused such that the motion of 
the system freezes out along these directions. The gas 
in the potential "dimple" provided by the laser then in- 
deed becomes an effectively one-dimensional condensate, 
in contact with the three-dimensional thermal cloud in 
the magnetic trap, which acts as its heat bath. The dy- 
namics of the order parameter is governed in this case 
by @J37| 
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' dt 



2y72 



2m 



$(z,i)+?7(z,t) 



where the external trapping potential in the weakly- 
confined direction V°^^{z) is again given in Eq. (63) and fi 
is the effective chemical potential of the one-dimensional 
system. The one-dimensional coupling constant g is re- 
lated to K by g — AiTKh^ /m. Physically, the function 
iR{z, t) describes the pumping of the one-dimensional 
condensate from the surrounding thermal cloud, and 
7y(z, t) corresponds to the associated noise with Gaussian 
correlations. Both these quantities depend on the one- 
dimensional Keldysh self-energy ?iE^(z), as discussed in 
detail in Ref. |^. For our purposes, we only need that 

iR{z,t) ^ -^nEK(z) 

2y72 



2m 



-f y'=''*(z)-M + T'"l^(2,t)rj , (89) 

{r^*{z, tUz\ t')) = ^S^(z)<5(z - z')5{t - t') , (90) 

where (...) denotes averaging over the realizations of the 
noise r]{z,t). The numerical techniques employed are 
discussed in Ref. where it also was shown that 

with the last two expressions, the trapped gas relaxes 
to the correct equilibrium, as ensured by the fluctuation- 
dissipation theorem. To simplify the numerics, the non- 
condensed part in the dimple is here allowed to relax to 
the "classical" value N{e) = [(3{e — and the com- 

parison to the previous mean-field predictions is there- 
fore carried out by making the same approximation in 
the calculation of both no(z) and n'{z). The normalized 
first-order correlation function at equal time g*-^^(0, z) 
corresponding to the previously computed phase correla- 
tion function, is calculated via numerical autocorrelation 
measurements, i.e., 



gW(0,z,i) 



($*(0,t)j>(z,t)) 



(91) 



where the brackets again denote averaging over the differ- 
ent realizations of the noise. Of course, the time t must 
be sufficiently large so that the gas has relaxed to thermal 
equilibrium and g^'^'^(0, z, t) is independent of time. 

In Figs. H and ^ we show the comparison of the many- 
body T-matrix theory to the above Langcvin calcula- 
tions, for the same temperatures used in Figs. ^ and 0. 
In Fig. I we compare the Langevin densities (|$(z,t)| ) 
to our classical mean-field density n{z). This yields ex- 
cellent agreement at low temperatures, except for a small 
region around the discontinuity in the mean-field theory, 
which can be understood from the fact that the local- 
density approximation always fails in a small region near 



the edge of the Thomas-Fermi radius. As expected, this 
region increases with increasing temperature. For T = 50 
nK, Fig. ^ further shows the deviation of the "classical" 
prediction of our mean-field theory from the "quantum" 



one calculated previously in Sec. VA and displayed in 
Fig. ^. Finally, Fig. |l^ shows the corresponding phase 
correlation functions as a function of position. Here we 
also find very good agreement in the entire temperature 
range. Note that the phase correlation functions are es- 
sentially indistinguishable for both classical and quantum 
treatments of the thermal cloud. 

It is interesting to note that the Langevin method 
yields continuous curves at the expense of computational 
time, due to the large number of independent runs that 
are required to reduce the statistical error. However, the 
Langevin method enables also a direct calculation of the 
time-dependent correlation properties via temporal auto- 
correlation measurements. Results of such studies, which 
are of the interest for the physics of an atom laser, will 
be presented in a separate pubhcation ]3^ ]. 

Finally, it is worth mentioning again that in obtain- 
ing our analytical expressions for the phase fluctuations 
and the density in sections II and V, we have used the 
many-body T-matrix for the interatomic interactions. As 



mentioned in section II B the many-body effects are im- 
portant in one and two dimensions. To appreciate this 
importance, we recalculate the density profiles and phase 
fluctuations using the two-body T-matrix. Thus for dis- 
tances below i?xF the differences are due to Eq. (pj]), 
whereas for distances above i?TF they are a result of 
Eq. (|2|). In Figs. ^ and |l| it is clearly seen that 
the inclusion of many-body effects has led to a better 
agreement with the exact Langevin results. Moreover, 
the many-body corrections become more pronounced at 
higher temperatures. In Fig. |l^, we show how the renor- 
malized interatomic interaction strength T^™{—2ii{z)) 
depends on position. We notice that the effects of this 
renormalization becomes most significant near the edge 
of the condensate and for temperatures closer to the 
transition temperature, as expected from the results of 
Refs. 
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VI. CONCLUSIONS 



The Popov theory suffers from infrared divergences in 
the equation of state at all temperatures in one dimen- 
sion and at any nonzero temperature in two dimensions. 
These infrared divergences can be traced to an inaccurate 
treatment of the phase fluctuations. We have proposed 
a new mean-field theory for dilute Bose gases, in which 
the phase fluctuations are treated exactly. We have also 
used this to arrive at an improved many-body T-matrix 
theory. The resulting equation of state is free of infrared 
divergences and the theory can thus be applied in any 
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dimension. Our modified many-body T-matrix theory 
is capable of reproducing exact results in one dimension 
and the results in three dimensions are to a large extent 
the same as those predicted by Popov theory. We have 
used the theory to calculate the reduction of the recom- 
bination rate of a spin-polarized two-dimensional hydro- 
gen system. Comparing our calculated rate with the ob- 
served values we found reasonable agreement, although 
more work is required to make a detailed comparison. 

We have also applied the theory to the Kosterlitz- 
Thouless phase transition. The modified many-body T- 
matrix theory is used to calculate initial conditions for 
the supcrfluid density and the fugacity of the vortices 
in a renormalization group calculation that incorporates 
the physics of vortex pairs. We have calculated the criti- 
cal temperature for a fixed value of the s-wave scattering 
length as a function of density, and it was found that 
Tc increases almost linearly with density. More precisely 
we have obtained nA^ ~ 7. Wc believe that this result 
gives a lower bound on the critical temperature, since 
the Kosterlitz-Thouless renormalization group equations 
do not include quantum effects, which in principle affect 
nonuniversal quantities. 

The modified many-body T-matrix theory was also ap- 
plied to calculate density profiles and phase correlation 
functions of a one-dimensional trapped Bose gas for a 
variety of temperatures. At very low temperatures, the 
phase correlation function was found to decrease only 
very slightly over the size of the system, indicating that 
the equilibrium state contains a true condensate. At 
larger temperatures, it decreases fastcir. and the gas now 
contains only a quasicondensate. In future work, we will 
look in detail at this and also at the full crossover prob- 
lem betwcxai one, two, and three dimensions. Finally, 
the densities and phase correlation functions predicted 
by our mean-field theory for various temperatures were 
compared to the corresponding predictions of a nonlinear 
Langevin field equation, which gives numerically exact 
results. The agreement was found to be very good for 
the entire temperature range studied. 
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FIGURES 
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FIG. 1. Reduction of the three-body recombination rate as a function of the density for a temperature of T = 190 mK and 
three different values of the scattering length. The dotted line corresponds to a = 2.4ao, the long-dashed line to a = 1.2ao, and 
the dashed line to a = 0.6ao, respectively. 



0.1 



0.08 



0.06 



0.04 



0.02 



0.25 



C=-0.0276 
C=-0.0278 




FIG. 2. Renormalization group flow for the coupling constants y and K. These curves are given by Eq. for different 
values of C. 
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FIG. 3. The critical temperature for the Kosterhtz-Thouless transition as a function of the density for spin-polarized atomic 
hydrogen with a = 2.4ao. 
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FIG. 4. The critical degeneracy parameter nA^ as a function of the density for spin-polarized atomic hydrogen with a = 2.4ao. 
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FIG. 5. Quasicondensate density no (solid curve) and superfluid density Us (long dashed curve) as a function of temperature. 
Also plotted are the Kosterlitz-Thouless condition nsA^ = 4 (dotted line) and the condition noA^ = 6.65 (dashed line). 
The Kosterlitz-Thouless transition takes place when the dashed line intersects the solid curve. At the intersection point the 
long-dashed curve reaches the dotted line. 




FIG. 6. Density profile of a trapped one-dimensional Bose gas at four different temperatures. The quantities h and k are 
defined in the text. 
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FIG. 7. The crossover temperature Tqc is shown with the sohd curve, and the chemical potential at this temperature is 
shown with the dashed curve, both as a function of the coupling constant. The temperature is scaled to To = /mk^ll. The 
inset shows the same curves on a double logarithmic scale. 




FIG. 8. Normalized first-order (phase) correlation function as a function of position for different temperatures. 
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FIG. 9. Comparison of the mean-field densities profiles (sofid curves) to numerical solutions of the Langevin equation in 
Eq. ds^ ) (noisy curves). All the above curves are calculated using the classical approximation of the Bose-Einstein distribution 
function. For the T = 50 nK case we have also plotted the corresponding density calculated using the full Bose-Einstein distri- 
bution function (dashed curve) in order to show the difference between the classical and quantum mean-field approximations. 
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FIG. 11. Study of the many-body renormalization effects on the density profiles. The exact results are also shown with the 
noisy curves. 
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